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Abstract 


It  Is   shown  that  the  differential  equation 

y'»     +      [X  +  $(x)]y     =     0 

can,  -under  suitable  conditions,  be  solved  by  ass\milng  a  solution 
of  the  form 

y     =     A(x)    sin  cp(x) 


where 


q)'(x)     =    \/rT"$(^+i^^^^|^  sin  2cp(x) 


^'W    =    -^WsTT^-^'^W 


Use  of  the  first  equation  leads,  when  boundajy  conditions  are 
applied,  to  asymptotic  estimates  of  the  eigenvalues. 

In  partictilar,  in  the  case  of  Hill's  equation,  it  is  shown 
that  the  Instability  intervals  vanish  faster  than  any  inverse 
power  of  k,  k  being  the  order  of  the  corresponding  eigenvalues, 
when  $(x)  is  an  analytic  function. 
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Introduction 
In  the  classical  SttLrm-Liouvllle  problem  one  Is  concerned 
with,  the  differential  equation 


^p(x)  g+  (Xp(x)  +  g(x))y  =  0 


subject  to  the  ■bo\xndary  conditions 

ay(0)   +    py'(O)      =     0 

7jM   +  5y'(rt)      =     0   . 

The  finite  interval  under  consideration  can  be  taken  to  be  the 
interval  (O^it) ,  X  is  a  parameter  and  p(x)  ,  p(x)^  g(x)  are  real 
functions,  and  furthermore  p(x)  and  p(x)  are  positive,   a,  p,  7,  5 
are  constants  such  that  |a|  +  |3|  ^  0  ^^^    l/l  +  |5|  4  0*   "-^^  ^^ 
then  interested  in  determining  those  values  of  X  for  which  non-trivial 
solutions  of  the  differential  equation  can  be  found,  which  satisfy 
the  boundary  conditions. 

If  p(x)  and  p(x)  are  doubly  differentiable,  one  can  by  means 
of  the  Llouville  transformation  reduce  the  equation  to  the  form 


,2 

^+  (X  +  $(x))y  =  0  . 
dx 

Borg'-  -J  discussed  this  equation  and  the  asymptotic  form  of 
the  characteristic  values  of  X.   His  chief  interest  was  a  discussion 
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of  the  inverse  Stiirm-Llouvllle  problem,  where  one  is  interested  in 
reconstructing  the  differential  equation  from  a  knowledge  of  the 
Sturm-Llouville  spectrum. 

In  this  publication  a  method  is  developed,  which  leads  to 
Borg's  results  in  a  very  efficient  fashion,  and  has  the  additional 
advantage  of  yielding  improved  error  estimates.   Furthermore  under 
fairly  general  conditions  one  can  easily  obtain  asymptotic  series  for 
large  eigenvalues  and  the  corresponding  eigenfunctions.   The  method 
leads  to  particularly  precise  information  in  the  case  of  Hill's 
equation. 

As  will  be  observed  subsequently  the  method  is  intimately 
related  to  the  W-K-B  method,  with  the  additional  advantage  of 
providing  a  rigorous  foiondation  for  the  latter  method. 

The  method 

If  we  consider  the  differential  equation 


y"  +  Q^(x)y  =  0 


where  Q(x)  Is  a  dlfferentlable  non-vanishing  function,  it  seems 
reasonable  to  assume  a  solution  of  the  form 


A(x)  sin  cp(x) 


If  in  addition  one  postulates  that 


y'  =  A(x)  Q(x)  cos  cp(x)  , 


-  5  - 

one  can  by  means  of  the  original  equation  show  by  a  direct 
calculation  that 


9 


'(x)   =  Q(x)  +  1^^  sin  2cp(x 


A'(x)   =   -A(x)  1^  cos^  cp(x)  . 

SimllsLTly,   if  one  assumes  that 

y     =     A(x)    cos  9(x) 

y'  =  -A(x)   Q(x)  sin  rp(x) 
one  obtains 

9'(x)   =  Q(x)  -  11^  sin  2  cp(x) 
A'(x)   =   -A(x)  l^  sin^rpCx)  . 

Evidently  once  'p(x)  has  been  found  from  the  first  of  these  equations, 
A(x}  can  be  obtained  by  an  integration  from  the  second  equation,  and 
it  also  follows  that  A(x)/A(0)  is  a  positive  fimction,  since  it  is 
an  exponential. 

The  Sturm-Li ouvi lie  problem 

We  now  turn  to  the  equation 

y"  +  [X  +  $(x)]y  =   0 

subject  to  the  boundaiy  conditions 
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ay(0)  +  Py'(O)   =  0 
7yM   +   5y'(3t)   =  0  , 
and  subdivide  the  problem  into  foiir  distinct  subcases: 

la.  a  =  7  =  1 

p  =  5  =  0 
Ila.  a  =  5  =  1 

P  =  0 

7  arbltraiy 

lb.  3=5  =  1 

a,  7  arbitrary 
lib.  p  =  7  =  1 

5  =  0 

a  arbitrary 

Problem  la. 


We  now  consider 


y"  +  [X  +  $(x)]y  =   0 


where     y(o)   =  y(jt)   =  0  , 

and  <J>(x)  is  a  differentiable  function  of  mean  value  zero^  i.e. 


$(x)dx  =  0 
o 
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Under  the  assumption  that 

y  =  A(x)  sin  cf(x) 


y'  =  A(x)  \/X  +   $(x)  cos  cp(x) 


we  obtain 


1  $'(x) 


cp'(x)   =   v^X~r¥(^  +  ^  ^^^^^1^^  sj,n'2rp(x 


where 

cp(0)   =  0 
9(rt)   =   fcrt 

and  k  is  any  sufficiently  large  integer.   Here  sufficiently  large 
means  so  large  that  for  the  corresponding  value  of  X 

\  +  $(x)   >  0  . 

The  differential  equation  can  be  rewritten  as  an  Integral  equation 

X  x 

9(x)   =  J    \/\  +   $(x)"  dx  +  i  J   ^  I  '^^^^  sin  2cp(x)  dx  . 
o  o 

From  this  equation  we  see  that 


^i)(x)   =   /  ^X  +  ii^   dx  +  O(-)  , 
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and  by  an  iteration 

X  X  ^   X 

cp(x)      =    J     VTTTT^  dx     +     I    J     Y^  sin  I2J   \/X  +  <I.(x)    dxj  dx 
o  0^0 


+      0(-^)    . 
X 

¥e  can  get  estimates  for  X  for  large  values  of  k  from  the  boundary 
condition 

Jt  rt  ^   x  ^^^ 

J     \/X  +   $(x)  dx  +  ^  j   ^  ^  '^^^^  ^^""[^f   ^^^^^^^  '^j  <^ 
o  o  ^  o  ^ 

+  O(^)  =  kjt  . 
X 


If  $'(x)  is  differentiable  an  integration  by  parts  can  be  applied  to 
the  second  integral  above  and  one  then  finds  that 

V/X  =  k  +  -i^  +  0(4)  , 


where 


31 

1  r.2,    .    ^  <!>'M    -  O'Co) 


or 


X  =  k^  +  ^  +  Oi±) 
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If,  on  the  other  hand,  <I>'(x)  is  not  dlfferentlable  ve  expand 
*(x)  In  a  Fourier  series  hy  continuing  il>(x)  outside  its  interval 
of  definition  (0,n)  as  an  even  function  of  period  2jt  and  mean  value 
zero.   Then 


$(x)  =   ^  a  cos  nx  , 


n 


and  we  find  that 


V^  =  k  +  ^  +  -^  f^^i^)   dx  +  0{-\) 


8itk' 

o 


Now  that  an  asymptotic  estimate  of  \  is  known  one  can  find  asymptotic 
estimates  of  cp(x)  and  A(x)  and  therefore  y(x) .  Then  one  finds  that 


y(x)  = 


X 

A(0)  sin  [2kx  +  ^  /  f  (x)  dx] 
o 

[k  +  $(x)]^/^ 


[x.o(i)] 


It  is  interesting  to  compare  the  eigenvalues  with  those  of 
the  corresponding  problem  for  $(x)  =  0.   Then 

y"  +  Xy  =  0 
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y(0)   =  y(rt)   =   0 


and 


\A  -  V^  =  A  +  o(-^) 


for  ^'(x)  dlfferentlable.   When  <i>'(x)  is  not  differentlable 


^^     8rtk^  -*  k^ 

o 


Problem  lb. 


For  the  equation  with  constant  coefficients  we  have  in  this 


case 


y "  +  y.y    =  0 
ay(0)  +  y'(0)   =  0 

7jM   +   y'(rt)   =  0 

so  that 

(7-a)  V^ 
tan  VX  «  = 


X  +  7a 
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One  can  show  from,  this  transcendental  equation  that 


c 

\ 

2 

^ 

^ 

7- 

a\ 

1 

7   -< 

X- 

7-a 

\  Jt 

-) 

+ 

3 

Jt 

irk 

1.5 

V5C  .  ^  .  ^  .  ^"-^    ^   '   .  o(-^) 

When  we  t\im  to  the  more  general  problem. 

y"  +    [x  +  ^(x)Jy     =     0 
ay(0)   +  y'(0)      =     0 
/(it)   +  y(it)      =     0 

we  assume  a  solution  of  the  foim 

y  =  A(x)  sin  (^(x)  +  A) 
y'  =  A(x)  Vx  +  $(x)  cos  (cp(x)  +  A)  . 

If 


_^  Vx  +  <Ko) 

A  =  -tan   


the  boundary  condition  at  x  =  0  is  satisfied.  As  before  we  find 


cp'(x)   =  MX  +   <5(x)   +  ^^^^^^  sin  [2cp+  2a] 
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vith.  the  "bouridary  conditions 


9(0)     =     0 


cp(rt)     =     krt  +  tan"       -     tan 

a  7 


As  In  Problem  la.  we  can  rewrite  this  equation  as  an  Integral 
equation  and  obtain  the  asymptotic  solution 

X  X 


sin    [2    /     \/X  +  <I>(x)    dx     +     2A]    dx     +     0(-^) 


Applying  the  boimdary  condition  at  x  =  Jt  we  find 


It  Jt  X 

J  VX  +  *(x)    ^  +  1    J   rT%^  sln[2  J    VX  +  *(x)    dx  +  2A]dx 
o  00 

1  _i  vx  +  $(0)       _i  v^rr^x^ 

+  0(-x)  =  kit  +  tan -  tan         r 


^2^ 


For  $'(x)  differentlable  we  can  solve  this  equation  asymptotically 
and  find 


-  u. 


^  =  -  -  s^ 


where 


'5    rt 


•-  o 


A  ccmparlsan  with,  the  corresponding  value  of  \/^   shows  that 


\/I  -  0( 


8k^      k^ 


When  *'(x)  is  not  differentiable  we  proceed  as  before  and 
expand  $(x)  in  a  Fourier  series.  We  then  obtain 


n 
/  0^(x)dx 


V^  -  v^ 


8k\ 


if  ^  °^^) 
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Problem  Ila. 

"She  same  methods  apply  to  this  piroblem-     The  results  merely 
vlll  therefore  "be  stated. 

The  equation 

y«'  +     Xy     =     0 

y(0)      =     0 

yy{n)  +  y'(n)     =     0 

has  the  eigenvalues  governed  by 


tan  V^  jt  =  


frcni  vhlch  one  can  find  that 


^  =  ^^^^-^^^-^^°<7' 


For  the  equation 


y'«  +  [x  +  *(x)]y  =  0 

y(o)  =  0 
yyM  +  y'(«)  =  0 
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we  have 


<p'(x)     =    0r-r^R    +    I    3-il^  sin  2cp(x) 


cp(0)     =     0 


<P(jt)     =     k  +  1  -  tan' 


^j_      \/X  +  *(rt) 


¥e  then  find,  for   *'(x)   dlfferentlable,  that 


0:  .    V^C     =       1     c^     +     0(4) 


8t 


3     2 


^? 


•where 


"     «  I  j      *^W<5x  -  M(rt)    -  *'(0)   -   <!"(«) 


When  *«(x)  is  not  differentlable 


Vx-  0C 


*  (x)dx 


8k' 


2k+l 


H^  +  f ) 


IT  -  0(-|) 


Prohlem  lib. 

Problem  lib  can  be  solved  Inrniedlately  by  realizing  that  under 
the  transformation 


li^ 


X  -*  rt  -  X 


7  -►  -a 


Problem  Ila  goes  over  into  lib  and  one  then  finds  for  ^'(x) 
different! able 


where  c,  =  —  < 
4   Jt 


$  (x)dx  +  1ki*(0)  +  'I''(0)  +  $'(jt) 


and  for  0'(x)  not  different! able 


\/l  -    \/k     = 


$    (x)dx 


8k^rt 


2k+l 


iv(kH-i) 


i:  ^  °^:i^ 


Hill's  equation 

We  now  consider  the  case  where  $(x)  is  an  even  differentiable 
function  of  period  rt.   For  boundary  conditions  we  use  the  following 
four  possible  sets  [see  reference  2] . 


1-)  y'(o) 

2.)  y(0) 

3.)  y'(o) 

h.)  y(o) 


y'(|) 


y(|)  = 


y(|)  = 


y'(f) 


=  0 
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which  lead  to  solutions  of  the  differential  equation  with  the 
following  character  respectively 

1.)  even  of  period  re 

2.)  odd  of  period  n 

5. )  even  of  period  2jt 

k.)  odd  of  period  2rt  . 

These  four  cases  correspond  to  Problems  Ib^  la,  lib,  Ila 
considered  in  the  previous  section,  respectively.   There  is  a 
slight  change  in  the  sense  that  the  interval  under  consideration 
now  is  the  interval  (O,  — )  ,  whereas  before  it  was  the  interval  (0,jt) 
The  only  change  that  must  be  made  in  cases  1.)  and  2.)  is  that 


9(|)   =  kJt 


We  then  find  for  cases  1.)  and  2.) 


VT  =  2k  +  — ^-   /  'I>^(x)dx  +  0(-4) 

o 


In  cases  5«)  aJi<i  ^O  we  have  to  modify  the  boundary  condition  at  jt 
to 


cp(|)   =   (2k  +  1)|  . 
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Then  we  find 


V^  =  2k  +  1  +    •*-,   /  $^(x)dx  +  O(^) 


o 


In  case  2.),   for  example,  we  can  show  easily  that 

cp(x)   =  2kx  +  \Jr(x) 

where  i|f(x)  Is  an  odd  function  of  period  :t.   If  we  use  this  form  of 
9(x)  as  a  trial  solution  in  the  equation 


1    1>'(x) 


9'(x)  =  \/X  +  $(x)   +  i  -^11^  sin  ap(x) 
cp(0)   =  0 


9(|)  =  k« 


we  obtain 

1    $'(x) 


i'(x)  =    \/x  +  $(x)   -  2k  +  ^  rr^fi)  ^^"^  t^^ ""  ^^""^^ 


t(0)  =  t(|)   =  0. 


From  this  equation  it  is  evident  that  \|/(x)  is  an  odd  function,  and  the 
equation  is  invariant  under  the  transformation 

x  -»  X  +  rt  . 
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Considering  the  boimdary  conditions  on  ^(x)^   the  result  is  evident. 

We  can  also  show  that  the   function  A(x)   has  period  n.      The 
differential  equation  for  A(x)    is 


*■«    =    -*('')    2(xl%))    -^'-pW 


S^     ■     -   h(xl''lU)      [l+   cos   (te+  2+W)] 


The  equation  is  invariant  \inder  the  transformation 


X  -♦  X  +  n 


Furthermore 


^""jm  ^  'if  X I  'ifx)  t  ^  ^°^  ^^  ^°^  2t(x) 

o 

-  sin  k'kx   sin  2i|f(x)   dx  , 

where  the  first  integral  vanishes  by  an  immediate  integration  and 
the  second  and  third  terms  vanish  because  they  represent  cosine  and 
sine  Fourier  coefficients  of  odd  and  even  functions.   Thus 


-  l[§}  =  0 
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so  that  A(jt)   =  A(0).   Thus  the  solution 

IT-    ■ 
y  =  A(x)  sin  cp(x)  .  .    , 

for  case  2.)  can  be  seen  to  be  aji  odd  function  of  period  jt  directly. 
It  is  easily  shown  that 


and 


tW   =  ^  j  $(x)  dx  +  0(-^) 


k 


so  that 


y  =  A(0)  sin 


1+kx  +  j^  j  $(x)  dx    +  O(^) 
o 


Evidently  these  asymptotic  estimates  can  be  Improved  If  necessary. 


The  stability  intervals  for  Hill's  equation 

It  is  veil  known'-  -^  that  if  the  eigenvalues  corresponding 
to  solutions  of  periodicity  jt  be  denoted  by  X   ,X    ,X^,    ...    and  those 
of  periodicity  2rt  be  denoted  by  A.',X',X',  ...  they  form  a  sequence 
ordered  as  follows 


^0  <  ^1  ^  ^2  <  ^1^  ^2'^H^K'^^^-\'^   • 


For  all  X  in  the  intervals(XQ,X^) ,  (X^,Xj)  ^    U^Ap  ,    (^^^ 
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the  differential  equation  has  only  bounded  solutions.   For  X  outside 
these  intervals  one  solution  is  unbounded.   Therefore  the  intervals 
(-"Aq)^  (X^A^).  i\Ar,),    {^,\),    •••  are  known  as  the  instability 
intervals  and  they  are  known  to  become  small  for  high  eigenvalues, 
but  little  is  known  regarding  their  order  of  smallness.   An  immediate 
consequence  of  the  asymptotic  formulas  developed  in  the  previous 
section  is  that 


k 


from  which  we  see  that 


1 


Vi  -  \  =   °b^  ' 


k^ 


and  similarly. 


K    =    o(-^) 


k+l   k     ^^y 


When  the  function  ^{x)    is  analytic  one  can  give  much  more 
precise  estimates  for  the  instability  intervals.   For  example,  for 
case  1. )  we  have 

:t/2  rt/2 

cp(|)   =  /  VX+  Mx)  ^-l    f      ^I'if^)    -i^  2cp(x)dx  =   k:t 
o  o 

and  for  case  2. ) 
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rt/2  rt/2 

9(f)  =  /  v^  +  *(x)  <ix  +  ^  /  X !  lu)  ^^"^  ^^""^  d^  =  fat 

o  o 

When  0(x)  is  analytic  so  is  il>(x).   Therefore  one  can  apply  Integration 

by  parts  to  the  integral 

jt/2  .       .■  It/2 

r  <I)'(x)         .      o,  ^    ^    ^  r  ^  ^    ^     (i  1  "^'(x)      ^ 

/        ^; 1/    \    sm  2cp(xj    dx     =      /    cos  2cp(x)   ^—  -r — pT^-y  T S    \    dx 

J        X  +  <|)(x)  ^^    '  J  ^^    ^    dx  2cp'(x)    X+  $(x) 

o  o 


It/2 


sin  2cp(x)  ^  ^^^.^  -  -^r^  rrifey  '^ 


:t/2 
J     ^^^  2cp(x)     -^j^,^  __^ax, 


where  the  integrated  terms  vanish  by  virtue  of  the  periodicity  properties 
of  the  integrand.   Since 

cp(x)   =  0(\/x)   =  0(k)   ''    ■  •■  ■     ' 

cp'(x)   =  O(V^)   =  0(k) 


it  follows  that 
«/2 
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for  all  values  of  N.  By  taking  the  difference  of  the  equations  for 
cases  1.)  and  2.)  we  see  that 


v^  -  v^  =  °  hk^ 


or 


\+l  ~  \    ~     °  TTlTl 


since  W  is  arbitrary  we  can  conclude  that  the  finstability  internals 

> 
vanish  faster  than  any  negative  power  of  k. 

A  similar  process  shows  that  when  $(x)  is  m  times  different! able 

^k+1  '  '^k  "  °  V^^l 

A  result  due  to  Erdfelyi. 

Erdelyi  derived  ^  -^    some  estimates  for  small  <I)(x).   These 
follow  very  simply  by  this  method.   We  write  e<Ii(x)  instead  of  $  (x) 
and  assimie  that  e  is  small.   Then  for  case  2.)  we  obtain 


cp'(x)   =  \/x+   e$(x)   +  ^  ~^-^||^  sin  2cp(x) 


cp(0)   =  0 

cp(|)   =  kjt  . 

We  now  try  to  find  a  perturbation  series  solution  for  cp(x)  of  the  form 
form      cp(x)   =   2kx  +   £>(;(x)   +   ••• 

and  for  X.  we  use 

X  =  hk^  +     en  +   .... 
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Retaining  only  €  terms  we  have  immediately 


^^      l6k 


Use  of  the  boundary  conditions 


i(0)   =  ^(f)   =  0 


and  the  Fourier  series 


0(x)   = 

>; 

n=l 

a  cos  2mc 
n 

shows  that 

M-  = 

^2k 
2 

so  that 

. v2     '^2k 
X  =  ilk   +  — 7:— - 


Similarly  for  case  1. )  we  find 


) v2     '^2k 
X   =   4k   -  — TT- 


so  that 
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l\+l  -^1   =   l^^2kl 


It  is  evident  from  the  previous  development  that  Erdelyi's 
result  is  not  correct  in  general,  hut  Magnus  has  shown'-  -'  that  in 
an  average  sense  it  is  correct.   If 


7^  <  /g  <  7^  < 


denote  the  \  corresponding  to  even  solutions  of  period  jt  and 


^  <  "2  <  ^3  < 


correspond  to  odd  solutions  of  period  «.   Magnus  shoved  that 


0     /    ,  ^  n  n       /  ,   2n 

n=l  n=l 


A  similar  result  holds  for  cases  3-  )  a-nd  k.  )  . 

A   Special  Case 

The  previous  results  were  all  derived  tmder  the  assumption 
that  '^(x)  was  differentiable.   One  can  show  that  the  asymptotic 
forms  developed  here  do  not  hold  when  <l>(x)  is  not  differentiable. 
Such  an  example  was  treated  in  detail  in  reference  [5] .   The  equation 
considered  there  is 

y"  +  X^Q(x)y  =  0 


21^ 


Q(x)   =   1 


X  <  1 


1  <   X   <  L 


and  it  Is  shown  there  that 


^2n 


2n 
1  +  a(L-l) 


2  ^  ^  2 


2n-l 


2n     ~[ 
1  +  a(L-l)J 


-  +  9  ^ 
2  ^  ^  2 


"■k 


2n-l    ~] 
+  a(L-l)J 


1  +  a(L-l)  i  2  ^  ^  2 


^2n-l 


2n-l    ']  £      £ 
1  +  a(L-l)   2  "^  ^  2 


[x]  denotes  the  greatest  integer  less  than  or  equal  to  x  and  the 
0's  must  satisfy  the  inequality 


0  <  0  <  1  . 


Here  evidently 


^1 


k  =  1  +  a(L-l)  *  °(^) 


and  the  instability  intervals  are  0(l), 
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